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, Configuration mixing calculations performed in terms of the Skyrme/Gogny Energy 

■ Density Functional (EDF) rely on extending the Single-Reference energy functional into 

^S) ' non-diagonal EDF kernels. The standard way to do so, based on an analogy with the 

pure Hamiltonian case and the use of the generalized Wick theorem, is responsible for the 

recently observed divergences and steps in Multi- Reference calculations. We summarize 

here the minimal solution to this problem recently proposed^ and applied with success to 
o 

particle number restoratiori^. Such a regularization method provides suitable corrections 
of pathologies for EDF depending on integer powers of the density. The specific case of 
fractional powers of the density^! is also discussed. 
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5t 1 1- Energy Density Functional methods 

The nuclear Energy Density Functional (EDF) method is a unique tool to study 
static and dynamical properties of nuclei in a unified framework^. Although the 
nuclear EDF shares several features with Density Functional TheorjEl (DFT), the 
strategy used is different as it embraces two successive levels of description. 

On the first level, traditionally called "self-consistent mean- field theory", 
Hartree-Fock (HF) or Hartree-Fock-Bogoliubov (HFB), a single product state $0 
provides the normal p'^" and anomalous k^" density matrices the many-body energy 
is a functional of. We call this method a single-reference (SR) EDF approach and 
denote by £sr[^o] = fsfl:[p°°, *] the actual EDF. Although such a restric- 

tion is not necessary, one usually builds the EDF from an effective vertex (of the 
Skyrme or Gogny type), whose parameters are adjusted to reproduce a selected set 
of experimental observations. Independently of the starting point, the EDF can be 
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2 Lacroix, Bender, Duguet 

written as in any arbitrary basis as 

i^kl ijkl 

,1 \^ „-,PPP rtOO„00 „00 , 1 -PKK nOOK-00*K-"" -I- 

' 6 / y ^ijklmn rli Hmj Hnk ^ 4 / , '^ijklmn ^li '^jk '^mn ^ ■ ■ ■ : 

ijklmn ijklmn 

where the first term accounts for the uncorrelated kinetic energy, whereas v'''', v'^'^, 
yPP'', . . .denote effective vertices associated with the different terms of the EDF. 
There are a few important comments to be made at this point. First, and although 
it formally resembles it, Eq. ([T]) should not be confused with the expectation value 
of a Hamiltonian containing two-body, three-body, . . . interactions in the Hartree- 
Fock-Bogolyubov state $o- For this to be true, specific properties, e.g. v^^f^i = v'^Jli 
and v!^j^.i = —v^^i). for all fc, I), would have to be satisfied, which is usually not 
the case in the EDF context. Second, most popular and performing EDF cannot 
be written under the form of Eq. ([1]) as they contain a dependence on a non- 
integer power of the (local) normal densit}H^. We anticipate, however, that future 
EDFs will be constructed under such a form, typically truncated at forth or fifth 
power. Indeed, the regularization procedure presented here is inapplicable to EDFs 
containing non-integer powers of the density matrice^. 

While static collective correlations, e.g. pairing and deformation, can be ac- 
counted for within the SR EDF formalism through the symmetry breaking of the 
auxiliary state $0i dynamical collective correlations requires to perform a so-called 
Multi- Reference (MR) calculations, traditionally denoted as "beyond-mean-field". 
Such an extension, built by analogy with the Generator Coordinate Method (GCM) 
in the Hamiltonian formalisnJ^, allows one not only to incorporate additional cor- 
relations but also to describe low-energy spectroscopy and transition probabilities 
between states characterized by symmetry-restored quantum numbers. In strict 
analogy with the Hamiltonian formalism, the MR EDF is written as 

. , ^ E{o.i}gMH. /o /i SmbI^o, $i] ($o|$i) 

^ ^" E{o,iKm/o/i(<i>o|<i>i) ' ^ ' 

where non-diagonal matrix elements (<I'o|i?|<&i)/($o|^i) have been replaced by 
their EDF counterpart £mr[^o,^i]- The weight functions / are determined either 
by symmetry considerations, by diagonalization, or both. The product states 
belonging to the MR set are chosen according to the collective modes one wants to 
describe. In the absence of a well-founded prescription to build £mr[^o, ^i], only 
specific constraints can be imposed. For a number of reasonJ^, it is necessary to 
impose that £mr['Pq, $o] = £sr['^o] and £mr[^i,^q] = £*mr{^o, Following the 
Hamiltonian formalism, the most natural guidance is provided by the generalized 
Wick theoremSj (GWT) which tell us that £mji[^q,^i] is obtained by replacing 
SR density matrices by transition ones, i.e. k"^, in Eq. ([1]). However, we 

have shown that the use of GWT-based functional energ y kern els is the source of 
the pathologies recently observed in MR-EDF calculation^^. 
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Configuration mixing within the energy density functional formalismipathologies and cures. 3 

1.1. Pathologies observed in configuration mixing calculations 

An example of deformation energy surface obtained through a MR calculation based 
on Particle-Number Restoration (PNR) is given for ^^O in Fig.[T]using the trilinear 
SIII Skyrme EDF. Starting from a SR-EDF built from the auxiliary state $o which 
explicitly breaks the particle-number symmetry, in order to account for static pair- 
ing correlations explicitly, dynamical pairing correlations associated with PNR can 
be incorporated through a MR EDF calculation. Building the MR set from product 
states rotated in gauge space 1$^) = e*'^*|$o), Eq. Q specified to PNR readP 

/ dip—^£MR[<^0,^^]{^0\%)- (3) 

In practice, Eq. ^ is numerically estimated using the FomenkJI^ discretization 
procedure. The left panel of Fig. [T] presents results obtained for two different num- 
bers of mesh points in the discretization (dotted and dashed lines) . Obvious patholo- 
gies are visible, i.e. (i) the estimate of the energy landscape does not converge and 
(ii) non-physical steps appear at particular deformations as one increases the num- 
ber of mesh points used in the Fomenko procedur^. Authors have not only also 
faced the problem for PNeEHIII 

but also when performing angular-momentum 
restoratiorJISI. It has been recognizecfm that divergences in PNR may appear when 
either a proton or neutron single-particle level crosses the Fermi en energy, as pair of 
states differing by tt/2 are orthogonal in this case, i.e. ($o|'I'ip) = 0. When the same 
(density-independent) vertices are used in the p-h and p-p channel and the exchange 
is properly taken into account, problems were shown to disappear. The pro blem has 
been characterized more precisely thanks to a complex plane analysisl^, demon- 
strating in particular the less obvious but more profound occurrence of steps. This 
technique, however, cannot be extended to the restoration of other symmetries and 
does not lead to a practical solution of the problem. 

2. Minimal solution to the problem 

We have recently shown that the origin of difficulties can be traced back to the 
strategy used to design energy kernels entering the MR-EDF, i.e. the use of the 
GWT as a guidance. An early hypothesis "'^"'^ has been confirmed latefl^by avoiding 
the use of the GWT in the PNR case. More recently, a general solution which 
applies to any type of configuration mixing has been formulatecffl. The technique 
makes use of the following trick: given a pair of quasi-particle vacua, denoted by 
|$o) and |$i) (with possibly ($o|<l'i) = 0), one can always find a simple "BCS like" 
expression connecting these two states, i.eW. 

= Coi n i^*pp + Kp ^t) ■ (4) 

'^The fact that no divergence occurs is due to the particular form of the functional used which is 
strictly bilinear in the same isospirffl 
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4 Lacroix, Bender, Duguet 

In the quasi-particle basis where Eq. (H]) is vahd, GWT-based energy kernels read, 
e.g. for a strictly bilinear EDF, as (omitting the kinetic term) 



where (/i, /l) denote a canonical pair in the specific quasi-particle representation, 
Zvu = {Bvv I A~^)* while and 4>^j, stand for the upper and lower components 
of the quasi-particle states^. Using expressions ([5][S]) is convenient as it separates 
the contributions remaining in the SR limit (line [5]) from the rest. It also allows 
to identify the different sources of problems, (i) Self-interaction: Well-kown from 
DFT, self-interaction relates to the fact that a particle should not interact with 
itself, which, however, happens when the vertices v^p are not anti-symmetrizecff^. 
This is almost always the case in actual EDFs at least because of the approxima- 
tions used to treat Coulomb exchange. Such a self-interaction, if present at the SR 
level, further contributes at the MR level (lines (HH?])). (ii) Self-pairing: this new 
conceplP is specific to EDFs treating pairing correlations explicitly and relates to 
the fact that two paired particles should not generate correlation energies, beyond 
their direct interaction, by scattering onto themselves. Again, such a spurious con- 
tribution appears at both SR and MR levels, (ill) Steps and divergences: As the 
energy kernel is multiplied by ($o|^i) c< A*^, in the MR energy (see Eq. ([2])), 
only terms with v ^ ^ ov v ~ p, va. line ([8]) can lead to divergences and steps when 
Al^ = 0. In the pure Hamiltonian case, i.e. ■uf^'J,; = v'^^j^i and v'^^^.i = -wf^'J^. for 
all {i,j,k,l), the dangerous contributions coming from the two terms in Eq. ([8]) 
exactly cancel out and no divergence or step occurs. However, when different or 
non-antisymmetrized vertices are used, as in the EDF context, divergences and/or 
steps are observed, as on the left panel of Fig. [T] 

The quasi-particle basis introduced above allows one to proceed to the analogy 
with the Hamiltonian formalism on the basis of the standard Wick theorem rather 
than on the generalized one. Comparing the results of the two schemes, one proveJ^^ 
that terms with v = ^ ov v = fi in line ([8]) should be zero in the first place and 
must be removed altogether. This not only solves problem (ill) entirely but also 
remove finite spurious contributions to the MR energy kernel. Such a regularization 
technique can be applied to any type of configuration mixing performed in terms 
of an EDF depending on integer powers of the densities. It has been successfully 
applied to PNbP, as is exemplified on the left panel of Fig. [1] using the SIII Skyrme 
EDF. The correction not only removes the dependence on the number of mesh points 
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Configuration mixing within the energy density functional formalismipathologies and cures. 5 




Fig. 1. Left: Uncorrected (dotted and dashed lines) and corrected (solid line) particle-number 
restored quadrupole deformation energy obtained for ^**0 with SIII and calculated with L = 5 and 
199 discretization points of the integral in gauge space. The two corrected curves are superimposed. 
Right: attempt to regularize the particle- number restored energy of obtained with SLy4 that 
contains a non-integer power of the (local) normal density. 



and the non-physical steps, but also corrects the energy landscape away from those 
steps. The case of EDFs depending on non-integer powers of the density matrix 
has also been analyze Although divergences can be removed using a variant of 
the method proposed in the complex-plane analysis demonstrates that the 

left-over fractional power p'' with < 7 < 1 is ill-defined as it generates cusps in the 
PNR energy landscape, (see right panel of Fig. [1]). Generally speaking, one cannot 
use a functional that is multi-valued over the complex plane. This has important 
consequences on the present and future of EDF methods. 
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